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Modelling the correlation between the thermal Sunyaev 
Zel'dovich effect and the cosmic infrared background 
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ABSTRACT 

We show how the correlation between the thermal Sunyaev Zel'dovich effect (tSZ) 
from galaxy clusters and dust emission from cosmic infrared background (CIB) sources 
can be calculated in a halo model framework. Using recent tSZ and CIB models, we 
find that the size of the tSZ x CIB cross-correlation is approximately 20 per cent at 
150 GHz. The contribution to the total angular power spectrum is around —2 /iK^ 
at £ = 3000, however, this value is uncertain by a factor of two to three, primarily 
because of CIB source modelling uncertainties. We expect the large uncertainty in 
this component to degrade upper limits on the kinematic Sunyaev Zel'dovich effect 
(kSZ), due to similarity in the frequency dependence of the tSZxCIB and kSZ across 
the frequency range probed by current Cosmic Microwave Background missions. We 
also find that the degree of tSZxCIB correlation is higher for mm x sub- mm spectra 
than mmxmm, because more of the sub-mm CIB originates at lower redshifts (z < 2), 
where most tSZ clusters are found. 

Key words: cosmic background radiation - cosmology: theory - galaxies: clusters: 
general - infrared: diffuse background - submillimetre: galaxies 
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1 INTRODUCTION 

The ability to connect the luminous baryonic matter our telescopes see with the underlying dark matter density field is 
fundamental to our understanding of large-scale structure and galaxy evolution. Combining the 'halo model' with a Halo 
Occupation Distribution (HOD) prescription achi eves this through assoc i ating luminous sources with collapsed dark matter 
haloes, the peaks of th e de nsity distrib ution (e.g., Peacock fc Smith 2000l : Scoccimarro et al. 2001 : Berlind &: Weinberg 2002 



Coorav fc Shethll2002l : also iBondll 19961 and references therein). This framework has been successfully applied to interpret the 



clustering properties of several source populations, including local galaxies from the Sloan Digital Sky Survey (SDSS), luminous 
red galaxies (LRGs), and high-redsh i ft Lyman-break galaxies, using mea s urements of the an gular correlation function (e.g., 



Zehavi et al. 2004; Ouchi et al 



20051 : 



Blake et al.ll2008l : IZheng et al 



20091 : 



Zehavi et al. 2011 



A natural application of the halo model is to understand the correlation between two differe nt source populations, each 
tracing the underlying dark matter but having different dependence on host halo properties (e.g., Coorav fc Sheth 20021 '). In 
this work we present such a treatment for the angular power spectrum arising from the c orrelation between unresolved clusters 
contributing to the thermal Sunyaev Zel'd ovich effect (tSZ; 



Sunyaev fc Zel'dovichlligTol ). and the dusty sources that make up 



the cosmic infrared background (CIB; e.g.. |Puget et al.lll996f ). Along with radio sources, these two popula tions are currently of 



inter e st as foregrounds in Cosmic Microwave Background (CMB) temperature anisotropy analysis (e.g.. lWhite fc Maiumdai 



2004 : iRighi et al. 



20081 ') . The tSZ effect arises from the Compton scattering of CMB photons by hot electrons in massive 
galaxy clusters, while the CIB sourc es are understood to be mostly high-redshift (z ~ 1 — 4) galaxies under going active, 



dust-enshrouded star formation (e.g.. Bond et al. 19861 . 199ll : Hughes et al. 19981 : Blain et al. 19991 : Drainel 20031 ^ . Absorpti 



of the starlight by the dust grains and subsequent thermal re-emission leads to strong rest-frame far-infrared emission, and 
cosmic expansion stretches these photons' wavelengths so they are detected in the mm-bands of contemporary CMB missions 
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including the South Pol e Telescope (SPT; 



Hall et all 



20icl) . the Atacama Cosmology Telescope (ACT: IPowler et all [ioiol ') , 



and the Planck satellite ( Planck Collaboration et al.ll2011al ). 

The size of the tSZxCIB cross-correlation provides physical information about the extent to which the CIB sources are 
associated with the clusters that c ontribute to the tSZ signal. Understanding how and why star form ation in massive systems 



I ^i-^i \r-2 I : II 1 

has evolved over cosmic time (e.g. IChapman et al.l 120051 : iPannella et al.l 120091 : iMagnelli et al.ll201in is key to understanding 



galaxy evolution. In addition. iReichardt et al.l (2012 - hereafter R12) have recently found that uncertainty in the tSZxCIB 
power significantly degrades constraints on the kinematic Sunyaev Ze l'dovich effect (kS Z) from SPT data. The kSZ arising 



20121 ). however improving constraints on 



from peculiar motion of galaxy clusters has recently been detected (|Hand et al. 
the kSZ power fr om patchy reionization - in which CMB photons scatter off electrons in expanding ionized gas bubbles 
( Knox et al. 19981 ) - is of great interest for constraining the duration of reionization and thus star formation in the early 



universe. For these reasons, a model for the tSZxCIB power is required, and we choose a halo model framework since, as 
stated, it provides a natural basis for relating two different tracers of the matter field, and, furthermore, it has already been 
used to study the tSZ and CIB anisotropics. 

We note that a similar halo model approach has r ecently been used to model the cross-corr elation of the tSZ effect with 
the integrated Sachs- Wolfe effect ( Taburet et al. 201ll ). and the distribution of galaxy clusters ( Fang et al. 20121 ). 

In Section 2 we collate the halo model equations for the tSZ and CIB power spectra and present expressions describing 
the tSZxCIB cross-correlation. In Sections 3 and 4 we combine recent models for the tSZ and CIB and calculate the tSZxCIB 
contribution as a function of frequency and angular scale. The uncertainties in current models and the role of upcoming data 
sets are discussed in Section 5 and conclusions follow in Section 6. 

For the results presented in Section 4 we adopt a flat, ACDM cos mology, with h = Hp /lOO km s~^ Mpc^^ = 0.70, 
fim = 0.27, fib = 0.045, and ag = 0.80, consistent with WMAP-7 results jKomatsu et aLlboill ). 



2 MODEL 

Analysis of mm CMB data is typically performed using thermodynamic temperature units, while in the sub-mm, where the 
CIB is dominant, it is more natural to work with units of flux density. The equations in this section contain a mix of units, 
however the plots in Section 4 are made in either flux or temperature units. The conversion from thermodynamic temperature 
to flux density requires multiplying by dBv{T) /OTIt^Tcmb ' where 

ae^, s_ 2fcBi^^V_ 

dT ^ ' c2 (e- - if ' ^ ' 

with X = hpiz/ksT and the present-day CMB temperature Tcmb = 2.725 K. Note that converting between units requires an 
integral across the bandpass filter for the finite bands of real experiments. 



2.1 CIB power 



The angular power spectrum of unresolved CIB sources at multipole moment £ from correlating two maps at frequency 
and 1/2 may be written as 



C 



CIB 



(2) 



where the clustered component, (^ciB-ciust^ arises from the correlation of images of different sources in the two maps, while the 
Poissonian shot noise component, C'^™^^, is due to the finite number of sources in any given observed field, and may be thought 
of as arising from the correlation of images of the same source in the two maps. For the range of angular scales considered 
(~arcminute and larger), t he CIB sourc es are effectively point sources, and so the CI B shot noise is scale-independent. 

In the small-sky limit (|Limbeilll95i ). the clustered power may be written as (e.g.. lKnox et al.ll200ll : lTegmark et al ] |2002l ): 



C 



d.,f(.) 



d/5 



z dz 

where x is the comoving distance, dVc/dz = x^dx/d2: the comoving volume element, d/*^'^/dz the redshift distribution of 
the CIB intensity (sometimes called dS/dz), and Pgai is the three-dimensional power spectrum of the dusty galaxies (and any 
other sources) that make up the CIB. 

We write Pgai as the sum of one- and two-halo terms, corresponding to the c ontribution from pairs of sources lying in the 
same dark matter halo, and different haloes, respectively ( Coorav fc Sheth 20021 ): 

Pgai (fc, ^) = P^ti {k,z) + P^ti {k,z), (4) 

where 



Pgti{k,z) 



dM^ — -2 Mgai(fc,M,z) 



'gal 



(5) 
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and 



P^ii(fc, z) = PoMik, z) ( / dA/ ^(M, z) h^{k, M, z) 



Ugal(fc, M, Z) 



(6) 



AM -''■■^'-'-'n^^.iz) 

The integral is over halo mass; drih/dM is the halo mass function, A'gai(A/, 2) the number of CIB sources hosted by a halo, 
«gai the Fourier transform of the spatial distribution of the sources within their halo, and 6h the halo bias relative to the dark 
matter power spectrum, Pdm, such that the halo and dark matter power spectra are related by 

P^{k,M,z) = bl{k,PI,z)PoM(k,z) 



The comoving CIB source number density, ?igai, is given by 

ngal(2) = 



dM^(Af,2){iVga,), 



(7) 
(8) 



and the angled brackets denote averaging over all haloes of mass M at redshift 2. 

Note that equation ((3| assumes that the spectral properties of the CIB sources are not dependent on properties of their 
host haloes. To make our approach more general, and to facilitate writing down the tSZxC IB cross-correlation term later, 
we a l low the CIB flux f rom a source at a given redshift to depend on halo mass (as in, e.g., Righi et al. 20081 : Sehgal et al. 



2OIOI : IShang et al.1l20l2l ) . This requires weighting each integral over halo mass with the mean source flux, {Sv{M, 2)), yielding 
(omitting dependence on angular scale, redshift and halo mass for brevity) 

/Ih 



^ {S.,){S.,){N,.l (iVgal-1)) ^2^^ 



dM 



and 



where 



p/2h 

- gal, 1^1 1^2 



= Pot. 



dnh (g.,)(jVgai) 

AMl -—— 6h %al 

dMi s^. 



... dnh , (S',2)(Afgai) 
dAf2 s„o 



S,{Z) = I dAf^(M,2)(S.(Af,^)>(iVga,). 

Substituting into equation (|3]), and recognizing that d/^™/d2 = s^dK/d2, gives 



with the following expressions for the one- and two-halo terms 

. dnh 



^CIB-2h I ^CIB-P 



^CIB-lh / , dVc 



d2 / ^^'^ ^ ^^"^^ (S.2>(A^gal (iVgal " 1)> ^al , 



and 



d2— — Pdm 

d2 



dAfl-^bh (5'.,)(iVgal)Ugal 



"'■"i"^ J Az'""'J dAfi 

The CIB shot noise contribution is given by ('e.g.. lScott fc White|[l999l ) 



dM2^&h (5.2)(iVga,) Wgal- 



^CIB-P 

so in terms of the halo model quantities deflned above we have 

^CIB-P 



AS 



CI 



Az^ f dAf^(S.,5.2>(iVgal>. 



(9) 

(10) 

(11) 
(12) 
(13) 
(14) 

(15) 
(16) 



Note that the shot noise depends on (S,y^Si/2) whereas the one-halo power depends on (S^i ) (5i/2 ) ■ If the framework described 
here is to be used to perform joint fits to the shot noise and clustered power, the distribution of CIB source flux at fixed M 
and 2 must therefore be considered. 

Further modifications to these equ ations, such as distinguishing between the first, 'central', source occupying a halo 
and subsequent, 'satellite', sources (e.g., Berlind et al. 2003 : Kravtsov et al. 2004 : Zheng et al. 20051 : Tinker et al. 2O10l ). may 
easily be made. 

The CIB intensity (one-point function), J^^™, is simply the integral of dJ5^^^/d2, and is therefore given by 



(17) 



Note that the clustered power, shot noise and CIB intensity may all have quite different redshift dependence, due to the 
different factors appearing in the integrals in these equations. This point is demonstrated in Section 4. 



2.2 tSZ power 

In the halo model, the tSZ angular power spectrum from correlating two maps at frequency ui and 1/2 is 



tsz 



^tSZ-lh , ^tSZ-2h 



(18) 
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We follow previous studies (e.g., Komatsu fc Selia^bood : Shaw et aLlbood : Battaglia et al.ir2O10l ) in treating each dark matter 
halo as hosting a single, extended source (electron population) for the purposes of contributing to the tSZ power. Unlike 
for the CIB sources, the one-halo term and the Poisson shot noise are then the same thing (which we will refer to as the 
one-halo term). The inter-cluster two-halo term has been ne glected in existing mm-band CMB analysis, since, due to the 
scarcity of massive clusters, it is important only for £ < 300 ( Komatsu fc Kitavama 19991 ). and on these scales the primary 
CMB anisotropy is completely dominant. 

Again invoking the small- sky approximation, we write the one- and two-halo tSZ components as (jKomatsu fc Kitavama 
19991 : iKomatsu fc Seliaklbood ) : 



C, 



-tSZ-lh 



and 



tSZ-2h 



dM^{M,z)ynM,z) 



dnh , 



dM -^{M, z) h^{k, M, z) y,(M, z) 



where Qu is the spectral function of the tSZ effect, given, neglecting relativistic corrections, by 

+ 1 



4 , with X 



fee 2c ^ 



and yi is the two-dimensional Fourier transform of the Compton Y-parameter, j/3d ( Komatsu fc Seliak 20021 ') : 



yi 



47rrs 



Axx^ yzu{x) 



sin(fe/4) . 



(19) 



(20) 



(21) 



(22) 



here x = r/rs, where r is radius from the centre of the (assumed spherically symmetric) halo, rs is a scale radius, and 
Is = Da/ts, with Da = x/(l + ^) the angular diameter distance to redshift z. 

Note that < for i/ < 217 GHz. This means that, for certain choices of and 1/2, the tSZ power is negative. This can 
also be the case for the tSZxCIB contribution. 

The integrated Compton Y-parameter from a single cluster (integrated over the projected area of the cluster on the sky) 
is given by 



Isz = 



dxx y3D{x). 



Summing the contribution to the total tSZ intensity from all clusters we then have 



/tS Z rj-i 



dz 



dK 

dz 



dM^.sz. 



(23) 



(24) 



2.3 tSZxCIB power 

We now apply the halo model to describe the tSZxCIB cross-correlation. There is a one-halo term arising from the correlation 
of CIB sources with the tSZ from their host halo, but there is also a two-halo term arising from the correlation between the 
CIB sources in one halo and the tSZ from another. Even if we consider the extreme limit in which there is no star formation 
and no CIB sources in those haloes massive enough to make a non-negligible contribution to the tSZ, this second term would 
still exist, provided there was some overlap between the redshift distribution of the CIB sources and tSZ haloes. This is a 
key prediction of this work. On large angular scales, where the two-halo term dominates, the tSZxCIB power is therefore 
relatively insensitive to the astrophysical processes that affect star formation in dense massive haloes. 
We write the total tSZxCIB power as 

^tSZxCIB _ ^tSZxCIB-lh . ^tSZxCIB-2h /r,r\ 

In terms of the quantities defined in Sections 2.1 and 2.2, the one- and two-halo terms are 

C£!l^l^2*^™""' = 2CMB j dz^ J dA'/^J/£^g^i(S'^2> +fl^2('S'^l)^(^gal>%al (26) 

and 

^tszxciB-2h ^ y^^^ I dz^PuM I dMi^&hy. I dM2^6h(p,,{5.2)+3.2(S.i))(iVgai)uga,. (27) 

The fact that, in the non-relativistic limit, the tSZ spectral dependence is not coupled to redshift or halo mass leads 
to these cross-terms being simpler than would in general be the case for two source populations. While equations 26 and 
27 include the tSZ and CIB contributions from both vi and V2, an interesting case to consider is that of two widely spaced 
frequencies, such that one of the two terms in parentheses in equations 26 and 27 is negligible due to the different frequency 
dependence of the tSZ and CIB. Constraints from cross-spectra at widely-spaced frequencies may actually provide the best 
opportunity for detecting and constraining the tSZxCIB term. This idea is discussed further in Sections 4.2 and 5.2. 
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Using this formalism, the tSZxCIB cross-correlation power requires no additional functions or parameters over those 
used to characterize the tSZxtSZ and CIBxCIB power individually. This is an advantage of the halo model approach over 
the template-based approach used to describe the tSZ and CIB power in R12, which requires the inclusion of a correlation 
coefficient ^ to specify the tSZxClB cross-power. Furthermore, the frequency and angular scale dependence of ^ may not be 
negligible, meaning fitting a single correlation coefficient may not be sufficient to describe the tSZxCIB power (see Section 
4.5). 



3 EXAMPLE CALCULATION 



Ifaving laid out an analytic halo model framework for the tSZxCIB power, we now use recent models for the CIB and tSZ to 
perform calculations at frequencies probed by current arcminute-resolution mm and sub-mm telescopes: ACT, SPT, Planck's 
High Frequency Instrument, and HerschePs Spectral and Photometric Imaging Receiver (SPIRE) instrum ent. We focus on 
150 GHz as the band primarily used to constrain small-sc ale CMB anisotropy power by ACT and SPT re.g.. lHall et~alll201(]| : 
Fowler et ailboiol : iKeisler et al.ll201ll : lDunklev et al ]|201ll ). and cross-correlations with higher frequencies as an opportunity to 
potentially improve constraints on the tSZxCIB amplitude and angular scale dependence. Specifically, we perform calculations 
using bandpass filters from the following frequency channels: SPT's 150 GHz (2100 /im), Planck's 857, 545, 353 and 217 GHz 
(350, 550, 840 and 1380 ^m), and SPIRE's 1200 GHz (250 /ini). Predictions are also made for the ACT 148 and Planck 
143 GHz channels. 



3.1 CIB model 



For the CIB clustering powe r we u se the model of Xia et al.l (2012 - hereafter X12). This work uses the CIB source model 
developed in Granato et al. (2004) and Lapi et al'l (2011 - hereafter Lll). Three populations of sources are considered: 
spiral galaxies, starbursts (these two populations contribute to the CIB for z < 2), and proto-spheroids (the dominant CIB 
sources at high redsh ift, which do not exist at z < 1). X12 perform a joint fi t to angular power sp ectra from Herschel /SPIRE 



sources at high redsh iit, which do not exist at z < 1). A12 periorm a joint h t to angular power sp ectra rrom Uerschel/bfltitj 
l Amblard etallboill) . Plan ck (jPlanck Collaboration et al.ll2011al ) and SPT jshirokoff et allboill ). using a simple HOD model 



based on lZheng et al.l (|2005l ). There are two free parameters in the model: a minimum halo mass for hosting a proto-spheroid 
source, A/min, and the power-law index, Osat, which describes the occupation of massive haloes with additional proto-spheroids. 
The results in Section 4 are shown for \og{Adn-iin /Mq) = 12.24 and asat = 1.81, the means of the marginalized parameter 
distributions obtained by X12. We find that the statistical uncertainty from varying these two parameters is not significant 
compared to other uncertainties in the tSZxCIB power (see Section 5.1). 

Note that the X12 model does not include any dependence of CIB source spectral properties (luminosity, SED, etc.) on 
host halo mass. Consequently, the angular power spectrum is evaluated simply using equation (3) and the d/^P'^/dz predictions 
from Lll. 



3.2 tSZ model 



For the t SZ we use the model des cribed in Efstathiou fc Migliacciol (2012 - hereafter EM12), which is based on the halo 
model of iKomatsu fc SeliakI \200i ). but with the empirical electron pressur e profiles from rece nt X-ray measurements of 
the Representative XMM-Newton Cluster Structure Survey cluster sample ( Arnaud et al. 2O10l ). The model has two free 
parameters: an overall amplitude, A (such that Cf'^ oc A), and an evolution parameter, e, which describes departures from 
self-similar evolution of the electron pressure profile with redshift. For the calculations in Section 4 we fix A = 0.95 and e = 0, 
which, for our choice of halo properties and bright cluster removal (see be low), correspon ds to £{£+ l)Cl^^ /2Tv\i^:iooo = 4 /iK^ 
at 150 GHz. This is consistent with recent ACT and SPT measurements (jPunklev et al.l 2011; R12). The effect of introducing 
departures from self-similar redshift-evolution in the electron pressure profile is discussed in Section 5.1.2. 



3.3 Halo definition 



We follow XI 2 in using the halo mass function from Sheth fc Tormen ( 19991 ) . and the associated large-scale halo bias from 



Sheth et al.l (|200ll ). This mass f unction approx i mates the number density of haloes as a function of halo virial mass, Mvir 



The electron pressure profiles of Arnaud et al. ( 2O10l ) and the EM12 model are expressed not in terms of Mvir, but M500, 
defined as the mass of the region within which the mean overdensity is 500 times the critical background density at that 
re dshift. When calculating ye we must therefore first calculate M500 as a function of Mvir and redshift (see, e.g.. Section 2.1 
of lKomatsu fc SeliakI [20021 ). 

Again following X12, we assume that the galaxy density profile, Mgai, equals that of the dark matter in an NFW halo 



( Navarro et al.lll996l ). udm, and that haloes are truncated at the virial radius, r^ii, for the purposes of hosting CIB sources. 
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As in EM12, however, contributions to the Compton Y-parameter are allowed out to 4rvii.. The effect of halo tru ncation 
radius on the tSZ and clustered CIB power spectra is not neglig ible (e.g.. lKomatsu fc Selia^bood : Iviero et al.lboogi '). and it 
is not clear that these assumptions are reasonable over the range of source populations and redshifts considered. In Section 
5.1.1 we consider the effect of forcing the CIB sources to lie at least some minimum distance from the centre of their host 
haloes, motivated by evidence that galaxies at the cen ters of groups and clu sters are typical ly un dergoing le ss-active star 
formation than th ose lying further out (e.g., [Kenmcutg 1983, iHashimoto et al.l 1998, iBai et al.l 2006, iBai et al.l 2007; see also 



Boselli fc Gavazzil 2006 and references therein) . 



3.4 Source removal 



Analyses of tSZ and CIB statistics are sensitive to the level of removal of bright resolved objects (clusters and point sources) 
from maps. Removal of objects down to the resolutio n limits of current instruments like Herschel /S PIKE, Planck, ACT and 
SPT has some effect on the tSZ one-halo power (e.g., Komatsu fc Kitavama 19991 : Shaw et al. 20091 ) and the CIB shot noise 
(e.g., Bethermin et al.l 2011 ), an d will also affect the clustered CIB power to a lesser extent (which is, as yet, not quantified; 



Addison et all 2012) 



see Section 3.1 of 

As stated, X12 use the Lll CIB source model, accounting for removal of bright CIB sources by ignoring contributions to 
dl^^^/dz from sources with flux above some threshold Scut, corresponding to the resolution limit of the relevant band (50 mjy 
for SPIRE bands, 6.4 mJy for the SPT 150 GHz band, and various values for Planck - see Table 3 of lPlanck Collaboration et al. 
2011a). X12 also, somewhat arbitrarily, discount contributions to the clustered CIB power spectrum from z < 0.25. We do 
likewise for the clustered CIB and tSZxCIB power; for the angular scales relevant for constraining the kSZ effect {I > 3000), 
the effect of this truncation is probably negligible, although this may not be the case on larger angular scales (see Section 
4.3). 

We remove the brightest 0.1 clusters per squa re degree in our model, roughly corresponding t o a level of cluster d etection 
currently achieved in deep ACT and SPT maps ( Staniszewski et al. 20091 : Marriage et al. 2011 : Benson et al. 2011 ). We do 
this to eliminate the effect of any massive, highly extended cluster at z ~ (which, should it exist, would surely be masked 
in any real power spectrum analysis). The exact level of cluster removal has little effect on the tSZxCIB power, especially for 
£ > 3000 (see Section 5.1.2). We define cluster brightness as the tSZ decrement integrated over the cluster's projected area. 
Due to effects such as the instrument beam and map filtering, this is not exactly the quantity that is measured in real maps, 
however, a more sophisticated treatment is beyond the scope of this work. 



4 RESULTS 

In this section we use the equations of Section 2 and the tSZ and CIB models described in Section 3 to calculate several 
quantities: the redshift-dependence of the tSZ and CIB intensity, the angular power spectrum of the tSZ, clustered CIB and 
tSZxCIB components for a range of frequencies, the redshift-dependence of the tSZ, clustered CIB and tSZxCIB power, and 
the tSZxCIB correlation coefficient ^ used by R12. 



4.1 Redshift-distribution of intensity, dJ/dz 

Figure 1 shows the redshift-distribution of the tSZ and CIB intensity, plotted as dlnZ/dz. We obtained dlnJ^^^/dz using 
equation (17). In these units the tSZ intensity is independent of frequency. While the tSZ effect integrated over the projected 
area of a cluster is independent of the cluster redshift (neglecting redshift-evolution in halo concentration), the scarcity of 
massive haloes at high redshift results in z < 2 clusters contributing the m ajority of the tSZ intensity. 

The CIB intensity redshift-distribution is shown for two recent models (|Bethermin et al.l 2011 - Bll, and lLapi et al.l 2011 
- Lll) and two frequencies. Increasing v and moving into the sub-mm results in an increase in the relative importance of 
low-redshift sources as well as an increase in the total CIB in t ensity (not shown in Figure 1). This beh aviour is common to 
other recent CIB source models (including Lagache et al. 20031 : Negrello et al. 2007 : Marsden et al.lEoilh . and arises from the 
thermal dust emission from CIB sources resembling a modified blackbody peaked at a rest-frame wavelength of ~ 100 ^m. In 
the Lll model, the peak in dJ'^'^/dz at z ~ 1.5 is due to starburst galaxies, whose contribution decreases rapidly at higher 
redshift such that for 2 > 2 the proto-spheroids are the sole contributors. The structure in the d/'^™/dz curves arises from 
summing the contributions from the spiral, starburst and proto-spheroid sources; the Lll 1200 GHz curve in Figure 1 is the 
sum of the three 250 /im curves plotted in Figure 3 of X12. 

The overlap between the tSZ and CIB curves on this plot give some indication of the size of the tSZxCIB cross-correlation, 
although, as shown in equations 26 and 27, the cross-power also depends on other factors, for instance the extent to which 
bright CIB sources occupy massive haloes for the one-halo term, and the effective redshift window function introduced by 
Pdm for the two-halo term. 
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Figure 1. Redshift-di s tribut ion of tSZ and CIB intensity. The tSZ intensity distribution is calculated from the model of 
lEfstathiou fc Migliaccicl 1I2OI2I) and, in these units, is in dependent of frequency. We show the CIB intensity distribution at 150 and 
1200 GHz predicted bv lBethermin et al.l (2011 - Bll) and lLapi et~all (2011 - Lll). The overlap between the tSZ and CIB curves gives an 
indication of the size of the tSZxCIB cross-correlation. The degree of correlation increases for higher CIB frequencies because dl'^^^ /dz 
becomes increasingly concentrated to low redshift. 



Consider cross-correlating a map at a fixed frequency where the tSZ makes a significant contribution (e.g., 150 GHz) with 
maps at higher frequencies. Figure 1 suggests that the degree of tSZxCIB cross-correlation will increase as the frequency of 
the second map increases (at least up to 1200 GHz). However, the degree of CIB x CIB correlation will decrease as d/*^^^/dz in 
the second map becomes more concentrated to low redshift. We therefore expect the importance of the tSZxCIB contribution 
relative to the cross-spectrum CIB power to increase as the separation of the two map frequencies increases. This is illustrated 
in Sectio n 4.2, below. 

The iBethermin et al.l (|201ll ) CIB model predicts far more low-redshift fiux than that of Lll. The results in Section 4 
follow X12 and use the Lll predictions, however we consider the effect of redistributing some of the high-redshift fiux to lower 
redshift on the tSZxCIB power in Section 5.1.1. 



4.2 tSZxCIB angular power spectrum 

Figure 2 shows the clustered CIB, tSZ and tSZxCIB contributions to the angular power spectra from cross-correlating a 
150 GHz map with maps at a range of higher frequencies. Note that the CIB Poisson shot noise is not shown because it is not 
yet well-constrained by data for the widely-spaced bands. Several features are apparent in these plots. Firstly, as expected, 
the size of the tSZxCIB power relative to the clustered CIB power generally increases as the frequency separation of the 
maps increases. The 150 x 353 and 150 x 150 GHz spectra are exceptions to this trend because both terms in parentheses in 
equations 26 and 27 are significant (in the other spectra, one of the two terms is highly sub-dominant). In the 150 x 353 GHz 
spectrum, the contribution from the tSZ in the 353 GHz map correlating with the CIB sources in the 150 GHz map acts to 
'fill in' around 10 per cent of the 150 tSZ x 353 CIB power, since the tSZ effect in the two maps has opposite signs. In the 
150 X 150 GHz spectrum, the significant presence of both tSZ and CIB fiuctuations in the map acts constructively and leads 
to a larger tSZxCIB component (relative to clustered CIB). 

All three components shown in Figure 2 consist of one- and two-halo terms. As stated earlier, the scarcity of massive 
clusters means the two-halo tSZ power is sub-dominant over the whole range of angular scales shown (although masking more 
bright clusters would increase the relative importance of the two-halo term). For the X12 CIB clustering model, the tSZxCIB 
one-halo power dominates the two-halo power at significantly larger scales than for the clustered CIB. This is mainly because; 

(i) the low-redshift spirals and starbursts responsible for two thirds of the tSZxCIB cross-correlation power (compared 
to a third or less of the CIB power) are hosted in low-mass, and hence low-bias, haloes (reducing the two-halo power - see 
equation 27), and 

(ii) there is no mechanism to limit the star formation or CIB source emissivity in massive haloes in the X12 model, possibly 
leading to an overestimation of the one-halo tSZ x CIB power. 
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Figure 2. Angular power spectrum from the tSZxCIB correlation predicted for cross-correlating a map at 150 GHz with maps at various 
other frequencies using the model described in Sections 2 and 3. The CIB power here is clustered CIB power only (the Poisson shot 
noise power is not shown). Note that the tSZxCIB cross-power component is negative for each case shown here, while the tSZ power 
is negative in each case except 150 X 150 GHz. The importance of the tSZxCIB contribution increases relative to the CIB contribution 
as the frequency separation increases because the CIB becomes more concentrated to low redshift at higher frequencies, increasing the 
overlap with the tSZ clusters (see Figure 1). The tSZxCIB component is uncertain by a factor of two to three; to avoid cluttering the 
figure we do not show uncertainty estimates here, but see discussion in Section 5.1. 



The relationship between CIB source emissivity and halo mass, and the modeUing of the CIB clustering, are the largest sources 
of uncertainty in our tSZxCIB power calculation (see Section 5.1.1). 



4.3 Redshift-distribution of power, dCt/dz 

Figure 3 shows the redshift-distribution of power, plotted as dlnC^/dz, for the tSZ, clustered CIB and tSZxCIB components. 
Comparing the left-hand panels with Figure 1, the tSZ power and intensity have similar redshift distributions, as do the 
clustered CIB power and CIB intensity at 1200 GHz. The distribution of clustered CIB power at 150 GHz (particularly the 
two-halo power) is concentrated to higher redshifts than the intensity, with only a small contribution from the peak in dl /dz 
at z < 1. In the X12 model, the proto-spheroid sources that contribute all the CIB at z > 2 are hosted in massive haloes, 
which at these redshifts are highly biased (increasing the importance of these sources for contributing to the two-halo CIB 
power but not the intensity). 

The redshift-distribution of the tSZxCIB power shows considerably less frequency dependence than that of the clustered 
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Figure 3. Redshift-distribution of the angular power spectrum, dlnCe/dz, for the tSZ and clustered CIB components (left panels; 
compare with dlnl/dz in Figure 1) and the tSZxCIB component (right panels). While the redshift-distribution of clustered CIB power 
changes significantly from 1200 to 150 GHz, that of the tSZxCIB power is fairly similar. The compl ex structure in d l n(7*^^'^^™/d2 i s 
due to the existence of three different populations of CIB sources in our adopted CIB source model jLapi et al.ll201ll) . Ixia et al.l hoi^ ) 
find that the clustering properties of the z < 1 CIB sources, which make a significant contribution to the tSZxCIB power, are not 
constrained by angular power spectra of CIB fluctuations from Herschel/SPIRE, Planck and SPT, meaning that the low-z tSZxCIB 
power is highly uncertain (see Section 5.1.1). 



CIB power. This is because most tSZ clusters lie aX z <2 where, according to the lLapi et al.l (|201ll ) model, the CIB sources 
are predominantly starbursts, and the shape of dl'^'^/dz is fairly similar for the range of frequenci es considere d here (Figure 
1). The peak at z ~ 0.15 in the right-hand panels of Figure 3 is due to local spiral galaxies in the lLapi et al.l (|201ll ) model. 
The sharp step at z = 1 occurs because the high-redshift proto-spheroid population comes into existence discontinuously; 
while unrealistic, this feature does not affect our analysis because the power spectrum remains smooth and the redistribution 
of sources required to achieve a gradual transition in ACi/dz will induce a small change in the tSZxCIB power compared to 
other modelling uncertainties (see Section 5). 



4.4 Comparison to primary CMB and other foregrounds 



Figure 4 shows the most recent 150 GHz ACT and SPT power spectra (|Das et al.l 2011: R12), with contributions from primary 
and secondary CMB temperature anisotropics and foregrounds overplotted. The tSZ x CIB power is a sub-dominant component 
but may be comparable to the kSZ at £ > 3000. Note that the kSZ power shown in the figure includes the contribution from 
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Figure 4. Comparison of powe r spectra of prim ary and secondary C MB temperature anisotropies and foregrounds at 150 GHz. The 
data points are the latest SPT llReichardt et all 2011; R12) and ACT l lDas et al.ll201l|) measurements; we simply overplot the various 
power sp ectrum com ponents here rather than performing a fit to these data. The CIB clustering power was reproduced from the 
model of IXia et al.l (20 11; X12), as described in Section 3.1. The tSZ power spectrum was obtained from the model described in 
lEfstathiou fc Migliacciol (2012; EM12), fixed to have l{l + l)Cf'^ /2tz\i^o,qqq = 4 /iK^ (see Section 3.2), and the tSZxCIB power, which 
is negative at 150 GHz, was calculated by combining the X12 and EM12 models, as described in Sections 2 and 3. We show the kSZ 
power calculated in ISehgal et al. ! 1I2OIOI '). Radio and CIB point source shot noise levels were taken from R12 and X12 respectively (the 
ACT data points have been corrected to account for the difference in radio source shot noise levels due to more sources being masked 
by SPT). The primar y lensed CMB power was calculated assuming a standard ACDM cosmology consistent with WMAP constraints 
jKomatsu et alj|201l[) . 



bulk electron motion in galaxy clusters and the intergalactic medium but assumes instantaneous reionization; including the 
effect of patchy reionization would increase this signal. Since the tSZxCIB power is negative for the principal CMB channels 
of ACT, SPT and Planck, we would expect uncertainty in the tSZxCIB power to degrade constraints on the upper limit of 
the kSZ. 

In principle, the tSZxCIB and kSZ components could be separated on the basis of their frequency dependence, however, 
we find that the frequency dependence is actually very similar across much of the frequency range probed by ACT and SPT. 
Figure 5 shows the frequency dependence of the tSZ, clustered CIB, tSZxCIB and kSZ power. The tSZ and clustered CIB 
power are - individually - easily distinguishable from a blackbody, however the tSZxCIB closely resembles a blackbody 
(horizontal line) for v < 200 GHz. This will further worsen kSZ constraints, and indeed R12 find that the kSZ upper limit 
is increased by more than a factor of two when the tSZxCIB correlation is allowed, despite using data from all three SPT 
channels. 

To assist in the analysis of small-scale CMB data, we have made the tSZxCIB curve from Figure 4 available to downloacf^. 



http://www.physics.ox.ac.uk/users/AddisonG/ 
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Figure 5. Frequency dependence of the tSZxCIB and kSZ power. The frequency dependence of the tSZ and clustered CIB power are 
shown for comparison; these signals are individually easy to distinguish from a blackbody but the tSZxCIB power has a frequency scaling 
that is very similar to that of the blackbody kSZ over the range of frequencies probed by ACT and SPT. 



4.5 Correlation coefficient 



While, as mentioned, an advantage of the halo model approach is that there is no need for introducing a cross-correlation 
parameter into our model, quantifying the size of the tSZ x CIB cross-correlation as a function of fre quency and angula r scale 
with such a parameter is worthwhile for the purposes of comparison with other models. R12 and IZahn et al.ll2012l (2012, 
hereafter Z12) define a correlation coefficient, , as 



c: 



/^tSZ ^CIB I 



(28) 



Note that Cf ]^^^ here includes both the Poisson and clustered CIB components. When any .^-d ependence i s negle cted, R12 



find (, = -0.18 ± 0.12 using SPT data a t 95, 150 and 220 G Hz. Z12 combine the tSZ model of Ishaw et al.l (|20ld l with the 
various CIB source models presented in IShang et al.l (|2012l l. and find —0.02 > ^3000 > —0.34, with decreasing on larger 
angular scales (see Figure 3 in Z12). 

Table 1 lists the values of at £ of 1000, 3000, 5000 and 9000 predicted from the X12 and EM12 models for cross- 
correlating a 150 GHz map with maps at the other frequencies shown in Figure 2. For the CIB Poisson shot noise levels, we 
follow X12 and use values of 6106, 4931, 1425, 175, 11 and 0.72 Jy^ sr~^ at 1200, 857, 545, 353, 217 and 150 GHz, respectively 
We do not attempt to quantify the uncertainty in ^£ (but see discussion of power spectrum uncertainties in Section 5.1). The 
~20 per cent correlation obtained for 150x150 and 150x217 GHz is consistent with the findings of R12, and suggests that 
^£ is not strongly dependent on frequency across the ACT and SPT bands. Physically this arises because the CIB is highly 



coherent at mm wavelengths, meaning that C^' 



and \/ C^™ in equation (28) have a very similar frequency scaling. 



It seems doubtful that a single correlation coefficient is sufficient to parametrize the tSZxCIB power for joint fitting to 
mmxmm and mmxsub-mm spectra. Table 1 shows an increase in of up to a factor of two for 150x1200 GHz compared 
to 150x150 GHz, which arises because the CIB redshift-distribution changes, becoming more concentrated to low-redshift, in 
the sub-mm bands (Figures 1 and 3). 

Including the Poisson, as well as clustered, CIB power in equation (28) is unsatisfactory since rare, bright CIB sources 
may make a significant contribution to the Poisson power but be almost wholly irrelevant for the tSZxCIB. This could be the 
case in the mm bands if, for instance, muc h of the Poisson p ower comes from the luminous proto-spheroids predominantly in 
M ~ 1012-5-13 M0 haloes at 2 > 1.5 ( as in lLapi et al.ll201ll). T able 2 shows our predictions for if we take C^™ in equation 
(28) to be the clustered CIB power only, as in lShirokoff et al.l (201l|). Unfortunately, due to the different ^-dependence of the 
tSZ, clustered CIB and tSZxCIB power, this leads to an increase in the ^-dependence of 

Our results suggest that the frequency and angular scale dependence of ^ may not be negligible, although further work 
is required to assess the extent to which fitting for a single correlation coefficient may bias other parameters. 
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Table 1. Degree of tSZxCIB cross-correlation, §f , for cross-correlating a 150 GHz map with maps at various frequencies 



t SPIRE Planck Planck Planck Planck 150 GHz'" 

1200 GHz 857 GHz 545 GHz 353 GHz 217 GHz 



1000 -0.27 -0.20 -0.14 -0.09 -0.11 -0.12 

3000 -0.34 -0.26 -0.20 -0.14 -0.15 -0.17 

5000 -0.35 -0.28 -0.22 -0.16 -0.16 -0.18 

9000 -0.33 -0.28 -0.24 -0.17 -0.17 -0.18 



"calculated with SPT 150 GHz bandpass filter; multiply values in this column by 1.1 or 1.15 for ACT 148 and Planck 143 GHz bands 

respectively 



Table 2. As Table 1 but with excluding the Poisson shot noise GIB power in equation (28) - considering only clustered GIB 



I SPIRE Planck Planck Planck Planck 150 GHz 

1200 GHz 857 GHz 545 GHz 353 GHz 217 GHz 



1000 -0.30 -0.22 -0.15 -0.10 -0.12 -0.13 

3000 -0.44 -0.36 -0.27 -0.18 -0.22 -0.23 

5000 -0.51 -0.43 -0.33 -0.23 -0.26 -0.26 

9000 -0.60 -0.52 -0.40 -0.28 -0.30 -0.29 



5 DISCUSSION 

5.1 Sources of uncertainty in tSZxCIB power 

Numerous sources of uncertainty in the tSZxCIB power are discussed below. Overall, we find that the contribution to the 
total tSZ X GIB amplitude uncertainty is around a factor of two from uncertainties in the GIB modelling, tens of per cent from 
the tSZ treatment, and 10 per cent or less from other uncertainties. 



5.1.1 CIB 



Two key assumptions in the X12 model are that (i) there is no connection between halo mass and GIB source spectral 
properties, and (ii) the clustering of GIB sources can be described by the following parameters: {Afmin, Afsat, (T(logio M),asat}, 
the parameter space that has been used to perform fits to angular correlation function measurements of (for instance) LRGs 
and SDSS galaxies (jZheng et allboosl : [Zehavi et al.lboilh . What kind of uncertainty could be introduced in our tSZxGIB 
calculation if these assumptions are incorrect? This is a difficult question to answer without re-fitting to the angular power 
spectra used by X12 with different assumptions and parameters. Here we attempt to quantify the uncertainty with several 
simple examples, focussing on the 150x150 GHz spectrum on angular scales of ^ ~ 3000, most relevant for attempts to 
constrain the kSZ power. 

Firstly, we note that 60 per cent of the tSZxGIB power comes from haloes of mass over lO^** A^q, compared to only 
10 per cent of the clustered GIB x GIB. If star formation and far-infrared dust emission in these massive haloes is strongly 
suppressed compared to in less dense environments then the one-halo tSZxGIB power could be lower by a factor of two with 
little impact on the clustered GIB power. 

Secondly, X12 find a^ax ~ 1.81 ± 0.04 for the proto-spheroid sources that dominate the clustere d GIB power spectra ; 



a preference for Qsat > 1 has als o been reported in previous (similar) GIB clustering analyses (e.g., lAmblard et al.l 12011 



Planck GoUaboration et al .l2011al ), corresponding to the number of sources occupying a massive halo somehow increasing faster 



than the mass of the halo. This relation is not ge nerally supported by simulations or analyses of clustering properties of other 
source populations, whic h favour ctsat ^ 1 (e.g., Kravtsov et al. 20041 ; Zheng et al. 2005 ; Reid fc Soergell 20091 ; White et al 



201ll ; IZehavi et al.ll2011f ). We consider the possibility that the HOD parameters used for the GIB analyses, while capable of 
yielding a good fit to the angular power spectra, do not actually describe the true source properties, and that in fact Osat — 1. 
Setting ftsat = 1 causes the clustered GIB power to increase by ~30 per cent, but the tSZxGIB power to decrease by 20 per 
cent, because reducing Qsat for fixed d7*^'^/d2 leads to more of the GIB contribution (both intensity and power) coming from 
low-mass haloes. This suggests that an inaccurate or incomplete HOD model may, potentially, be causing an overestimate of 
the tSZxGIB at the tens of percent level. 

The results in Section 4 assume that GIB sources trace the dark matter density distribution in their host haloes within 
one virial radius (with the first, 'central', source sitting at the very centre of the halo), as in X12. If we force all the GIB 
sources to lie at radii r satisfying 0.5 < r/rvir < 1.0 (though still tracing the dark matter), the clustered GIB x GIB power at 
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£ = 3000 decreases by only a few percent compared to 16 per cent for the tSZxCIB power. With a more extreme exclusion 
of CIB sources from the central parts of haloes, such that 1.0 < r/vvir < 2.0, the clustered CIBxCIB power decreases by 15 
per cent and the tSZxCIB power by almost 60 per cent. 

As shown in Figure 1, there is significant scatter in the low-redshift CIB flux predicted by current CIB source models. 
We find that increasing the Lll low-redshift spiral and starburst flux by 50 per cent (with a rescaling of the proto-spheroid 
flux such that the 150 GHz clustered CIB power at £ = 3000 is unchanged) leads to a 30 per cent increase in the tSZxCIB 
power. If the low-redshift source flux contribution is instead doubled (again with a rescaling of the proto-spheroid flux), the 
tSZxCIB power increases by some 60 per cent. 

We flnally note that X12 treat the spiral and starburst sources (SS) and the proto-spheroid (PS) sources completely 
independently. In fact, since both populations exist for 1 < 2: < 2 (see Figure 3 of X12), with a significant overlap mass- wise 
in their host haloes, there should also be SSxPS contributions to both the one- and two-halo clustered CIB power, which X12 
do not include in their fitting. 



5.1.2 tSZ modelling and data 



The amplitude of the measured tSZ power spectru m is currently uncert ain at the 10-30 per cent level, depending on what 
is assumed about, for instance, the kSZ component (jPunklev et al.ll201ll . R12). This corresponds to a smaller uncertainty in 
the amplitude of the tSZxCIB power, which is subdominant to the CIB modelling uncertainty discussed above. 

As stated, the para meter e in the EM12 model describes departure from self-similar redshift evolution of the cluster 
electron pressure profile. Planck Collaboration et al. ( 2011bh find e = 0.66 ± 0.52 for a combined sample of clusters from 
Planck and the Meta-Catalogue of X-ray detected Clusters (MCXC) at z < 1 (see their equation 8 and Table 6, and also 
Section 2 of EM12). Holding the amplitude of the tSZ power aX I = 3000 fixed, we find a reduction in C^igooo™ of and 
30 per cent for e = 0.5 and 1.0 respectively, compared to the e = case. Note that holding C^£,|ooo fixed requires setting the 
amplitude parameter A = 1.24 and 1.58 for e = 0.5 and 1.0 respectively. This redshift-evolution uncertainty is probably the 
largest source of uncertainty in our tSZxCIB model after the CIB modelling discussed above. 

The results in Section 4 are shown with the brightest 0.1 deg~'^ clusters removed. For our adopted cosmology and halo 
definitions, this level of removal mainly affects pow er on large angula r scales, with a negligible effect in the tSZ and tSZxCIB 
power at ^ > 3000 (consistent with the findings of Shaw et al. 20091 ). We investigate the effect of removing ten times more 
clusters (1.0 deg'^) and find that, while the tSZxtSZ power is reduced by almost a factor of two at ^ = 3000, C^^zxCib 
decreases by only ~10 per cent. As might be expected, removing more bright clusters increases the size of the tSZxCIB power 
compared to the tSZ. 



5.1.3 Other uncertainties 



(i) Relativistic correct i ons a nd bulk motion of cluster electrons are expected to modify the tSZ power by ~10 per cent 
(e.g., Nozawa et al. 19981 . 20061 ). and will have a reduced effect on the tSZxCIB power. 



(ii) Statistical uncertainty in proto-spheroid clustering parameters from X12 contributes a 10 per cent uncertainty at 
£ = 3000 in the clustered CIB power (inducing an uncertainty of ~5 per cent in the tSZxCIB power), mainly due to 
the measured power spectra becoming dominated by shot noise on small scales. This measurement uncertainty is currently 
significantly smaller than the CIB modelling uncertainty discussed above. 



(iii) Uncertainty in ACDM cosmological parameters: the tSZ effect depends very strongly on (Tg - C, 



tSZ ^ 



=300 



model (C^^,3QQQ oc cr^ ); we find that changing (Tg by 0.03 (corresponding to the Icr WMAP-7 uncertainty; 



PC gs lu our 



Komatsu et al. 



20 111 ) leads to a ~20 per cent change in the tSZxCIB amplitude, however, if the tSZ amplitude is held fixed, C^^3qq*jj'^ changes 



by < 10 per cent. If a standard cosmology is assumed, uncertainty in cosmological parameters is not currently a dominant 
source of tSZxCIB uncertainty. 

(iv) We have not investigated varying halo properties (e.g., mass function, concentration, or introducing scale-dependent or 
stochastic halo bias). We postpone assessing the impact of uncertainties in these quantities on the tSZxCIB power to future 
work given the size of the CIB modelling uncertainties discussed above. 



5.2 Future data 

5.2.1 Angular power spectra from mm and suh-mm maps 

Our results show that the contribution of the tSZxCIB power relative to the clustered CIB plus tSZ power increases for 
cross-correlating mm and sub- mm maps with widely-spaced frequencies (Figure 2). Cross-correlating ACT and Balloon-borne 
Large- Aperture Submillimeter Telescope (BLAST; observations at 600, 860 and 1200 GHz) ma ps of an 8.6 deg^ a rea of sky 



near the south ecliptic pole revealed a significant detection of cross-spectra clustered CIB power jHaiian et al.ll2012l '). however 
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the small sky area means these data are unable to provide constraints on the tSZxCIB power. Imminent Herschel /SVYKE 
observations, also at 60 0, 860 and 1200 GHz, will provide ~100_deg^_of overla p with existing AClUI and SPlQ fields (as well 



as upcoming ACTPol - iNiemack et al.l 2010 - and SPTpol - iMcMahon et al.l 2009 - fields), which should allow us to improve 



constraints on the GIB and tSZxCIB correlations. 

In Section 4 we showed results calculated using the SPT 150 GHz bandpass filter. Since the GIB intensity and clus- 



tered power amplitude are falling off very steeply with decreasing frequency in the mm-bands (e.g. Fixsen et al. 19981 : 



Planck Collaboration et al.l l2011al l , t he contribution to t he power spectrum from clustered GIB sources is 
lower in the AGT 148 GHz channel (|Addison et al.ll2012l ) 



20 per cent 

This means that the tSZxGIB power in the AGT 148 GHz x 
SPIRE cross-spectra will be around 10 per cent larger relative to the GIB x GIB contribution compared to in SPT 150 GHz x 
SPIRE. The SPT 95 GHz channel is also potentially well-suited to constraining the tSZxGIB signal; the ratio of the tSZxGIB 
to the GIB power in SPT 95 GHz x SPIRE spectra wiU be some fi ve times larger than for SPT 150 GHz x SPIRE, based on 
the frequency scaling of the tSZ and clustered GIB (jAddison et al.l 2012; R12). Gross-correlations of mm and sub-mm Planck 
High- and Low- Frequency Instrument channels (e.g., 70x857 GHz, 100x857 GHz) may also be of use for constraining the 
tSZxGIB power for I < 2500. Which channels give the best constraints in practice will also depend on the amount of noise in 
the maps and the size of the Poisson GIB power. 



5.2.2 Other statistics 

The angular power spectrum of mm and sub-mm maps is not the only statistic with the ability to explore the connection 
between clusters and GIB sources. Other measurements that may also provide insight include: 

(i) cross-correlations of GIB-dominated maps with X-ray cluster maps (either stacking analyses or power spectrum based), 

(ii) targeted imaging of clusters in sub-mm bands (using, e.g., Herschel, Submillimetre Gommon-User Bolometer Array-2 
and the Atacama Large Millimeter/sub- millimeter Array) - this will improve constraints on the spatial distribution of GIB 
sources within massive haloes, which enters into our modelling through Wgai, and 

(iii) cross-correlations of GIB maps with catalogues of objects at well-measured redshifts (e.g., LRGs, quasars) - this will 
improve constraints on d/^^'^/dz, particularly at z < 1, most relevant for the tSZxGIB power; AGT, Herschel, AGTPol and 
Planck (as well as, e.g., IRAS at higher frequencies) have - or will have - the overlap with optical surveys required for these 
studies. 



6 CONCLUSIONS 

We have presented equations describing the contribution to the angular power spectrum from the correlation of tSZ clusters 
and GIB sources in a halo model framework. We then used these equations to calculate the tSZxGIB power spectrum at 
150 GHz and for a range of cross-spectra using recent tSZ and GIB halo models. We find that: 

(i) The tSZxGIB is a sub-dominant component of the angular power spectrum at 150 GHz, contributing approximately 
—2 at £ = 3000; uncertainty in the tSZxGIB power will degrade kSZ constraints, as found in recent SPT analysis (R12, 
Z12), due to the similarity in their frequency dependence. 

(ii) The size of the tSZxGIB power relative to the clustered GIB power increases if we correlate mm and sub-mm maps 
with increasing frequency separation (at least up to 1200 GHz, the highest frequency _ffersc/ie//SPIRE channel). 

(iii) Uncertainty in the amplitude of the tSZxGIB power is currently a factor of two or three, with the uncertainty 
dominated by uncertainty in modelling of the clustering of GIB sources (although there are other significant sources of 
uncertainty, including possible evolution in cluster electron pressure profiles with redshift). 

The framework laid out in Section 2 may be used to perform joint fitting to AGT, SPT, Herschel, Planck and other data 
sets. The halo model and HOD formalisms are not without their complications, and work remains in order to find the best 
way to parameterize the connection between GIB sources and their haloes; it seems essential that a thorough exploration of 
modelling assumptions accompany future analysis. 

Detection of the tSZxGIB power in mm x sub-mm cross-spectra may enable us to constrain the tSZ effect from unresolved 
clusters in spectra that are free from the signal of the primary GMB (although the CMB fluctuations will still exist as noise 
in the mm maps). Whether this can lead to improved tSZ constraints over mm-band data alone will depend on the true size 
of the tSZ X GIB signal and future improvements in the modelling of the GIB sources. 



2 http://horschcl.esac. esa.int/Docs/A02/GT2_accepted.html#GT2_mvicro_l 

3 http://hcrschcl.csac.esa.int/Docs/A02/OT2_accepted.html#OT2_mvicro_2 
* http: / /herschel. csac.esa.int/Docs/AOl/OTl_accepted.html#OTl_jcarls01_3 
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We finally remark that the tSZxCIB power may be a significant contaminant for detection of signals beyond the kSZ, 
for example, detection of the Integra ted Sachs- Wolfe effect by cross-correlating mm-band CMB maps with tSZ or CIB maps 



( Taburet et al.ll201ll : lllic et al.ll2011 
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